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Abstract 

An adaptive isogeometric method based on d-variate hierarchical spline con¬ 
structions can be derived by considering a refine module that preserves a 
certain class of admissibility between two consecutive steps of the adaptive 
loop [6]. 

In this paper we provide a complexity estimate, i.e., an estimate on how 
the number of mesh elements grows with respect to the number of elements 
that are marked for refinement by the adaptive strategy. Our estimate is in 
the line of the similar ones proved in the finite element context, [3 Ell- 
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1 Introduction 

Throughout the last years, Isogeometric Methods have gained widespread interest 
and are a very active field of research, [ME], investigating a wide range of appli¬ 
cations and theoretical questions. Due to the tensor-product structure of splines, 
there exists very stable procedure to perform mesh refinement and degree raising 
which are known in the literature as h-refinement, p-refinement, /c-refinement HD]. 
While these algorithms are very efficient, the preservation of the tensor-product 
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structure at least locally to each patch, produces a dramatic increase of degrees 
of freedom together with elongated elements. Mainly for this reason, several ap¬ 
proaches have been proposed to alleviate these constraints and they all need the 
dehnition of B-splines over non-tensor-product meshes. Indeed, there are several 
strategies and we mention here T-splines pQ, hierarchical B-splines [T^lfT^I^ and 
THB-splines [IB], but also LR splines [121 E], hierarchical T-splines [2], modihed 
T-splines ra. PHT-splines [TTl [26] amongst others. 

Clearly, the development of adaptive strategies exploiting the potential of non- 
tensor-product splines is an interesting and important step which has been ap¬ 
proached in a number of papers, at least from the practical point of view. In 
fact, despite their performance in experiments [H [T3112112D1 E], the advantages 
of mesh-adaptive Isogeometric Methods have not been assessed in theory until 
today. Partial results, in particular on approximation, efficient and reliable error 
estimates, and convergence of the adaptive procedure, have been proven in pre¬ 
liminary work [B] in the context of (truncated) hierarchical splines, and we aim at 
continuing this study providing further ingredients that are needed to prove the 
optimal convergence of the proposed adaptive approach in the spirit of Adaptive 
Finite Element Methods 

In particular, in this paper we address the complexity of the mesh rehnement 
procedure proposed in [B|. The relation between the set of marked elements and the 
overall number of rehned elements introduced by the rehne module is not straight¬ 
forward: additional elements may be rehned to create only (strictly) admissible 
meshes. Admissibility is a restriction to suitably graded meshes that allows for 
the adaptivity analysis of hierarchical isogeometric methods. By starting from an 
initial mesh conhguration, the complexity estimate provides a bound for the ratio 
of the newly inserted elements and the cumulative number of elements marked 
for rehnement in the subdivision process that leads from the initial to the hnal 
mesh. This allows to control the propagation of the rehnement beyond the set 
of elements initially selected by the marking strategy. An analogous complexity 
analysis is currently available for bivariate and trivariate T-splines [221 |2T] . 

This paper is organized as follows. In Section |21 we recall notation and ba¬ 
sic results from [Bj. Section [3] is devoted to the announced complexity estimate. 
Conclusions and an outlook to future work are given in Section 01 

2 Hierarchical refinement 

In this section, we recall notation and basic results from [B]. Since the complex¬ 
ity analysis of the REFINE module can be directly performed in the parametric 
setting, we avoid to introduce the two diherent notations for parametric/physical 
domains (corresponding to with/without the hat symbol^in [B]). 
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2.1 The truncated hierarchical basis 

Let C C ■ ■ ■ C be a nested sequence of tensor-product d-variate 

spline spaces defined on a closed hypercube D in For each level with 
^ = — 1, we denote by the normalized tensor-product B-spline basis 

of the spline space dehned on the grid and we assume that G° consists of 
open hypercubes with side length 1. The Cartesian product of d open intervals 
between adjacent (and non-coincident) grid values dehnes a quadrilateral element 
Q of G^. For all Q E we denote by hq := 2“^ the length of its side, and by 
i{Q) its level, i.e., i{Q) = k. Moreover, we assume a hxed degree p = (pi,... ,Pd) 
at any hierarchical level. The analysis could be generalized to the more general 
case of a non-uniform initial knot conhguration by suitably taking into account 
the corresponding maximum local mesh size and to variable degrees as well. 

In order to dehne hierarchical spline spaces, we additionally consider a nested 
sequence of domains 3 3 3 that are closed subsets of D. Any 

is the union of the closure of elements that belong to the tensor-product grid of 
the previous level. The hierarchical mesh Q is dehned as 

Q := {Q€g‘,l! = 0,...,N-l], ( 1 ) 

where 

:= {Q E G^ : Q C Q* E G^\ £* > i : Q* C A Q* C Q} (2) 

is the set of active elements of level i. We say that Q* is a rehnement of Q, and 
denote Q* ^ Q, when Q* is obtained from Q by splitting some of its elements 
via “g-adic” rehnement. Although the hierarchical approach allows us to consider 
any integer g > 2, we will focus on the case of standard dyadic rehnement with 
g = 2. Hierarchical B-splines are constructed according to a selection of active 
basis functions at diherent levels of detail, see also USES]. 

Definition 1. The hierarchical B-spline (HB-spline) basis Ti with respect to the 
mesh Q is defined as 

n{Q) := [fiEB^ : supp /? C A supp fi ^ ^ = 0,..., - l} , 

where supp denotes the intersection of the support of with 

The following dehnition introduces the truncation mechanism, the key concept 
used to dehne the truncated basis for hierarchical splines [16]. 

Definition 2. Let 

s= ^ cJ+^(s)/3, 


3 


be the representation of s E C. with respect to the finer basis . The 
truncation of s with respect to is defined as 

trunc^’*'^ s := ^ c^'^^(s)/3. 

/3eB^+i 

supp 

Definition 3. The truncated hierarchical B-spline (THB-spline) basis T with re¬ 
spect to the mesh Q is defined as 

T{Q) := {Trunc^+V:/5eS'n?^(Q), £ = 

where Trunc^'*'^ fi := trunc^“^(trunc^“^(... (trunc^+^(/3)) ...)); for any fi E B^ D 

n{Q). 

For details on the properties of the truncated basis, we refer to [IS1[I3. 

2.2 Admissible meshes 

We consider the class of admissible meshes introduced in [6]. 

Definition 4. A mesh Q is admissible of class m if the truncated basis functions 
in T{Q) which take non-zero values over any element Q E Q belong to at most m 
successive levels. 

Since the case m = 1 refers to uniform meshes, we will from now on focus on 
the case m > 2. The relevance of admissible mesh conhgurations relies on two 
properties of THB-splines. First, for each element Q of an admissible mesh, the 
number of truncated basis functions of degree p = (pi,... ,pfi) which are non-zero 
on Q is less than + !)• Second, if Q is an admissible mesh of class 

m, then for all truncated basis functions r G T{Q) and elements Q E Q with 
Q n suppr 7^ 0, we have |(5| < |suppr| < |(5|, where the hidden constants in 
these inequalities depend on m but not on r, Q and N. These properties may be 
suitably exploited in the analysis of adaptive isogeometric methods. 

In order to characterize a certain class of admissible meshes, we consider the 
generalization of the support extension usually considered in the tensor-product 
B-spline case to hierarchical conhgurations. 

Definition 5. The support extension S{Q,k) of an element Q E with respect 
to level k, with 0 < k < i, is defined as 

S(Q, k) := {Q' E G^ E B\ supp fi H Q'^ 0 A supp fi n Q ^ 0} . 
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By a slight abuse of notation, we will also denote by S{Q, k) the region occupied 
by the closure of elements in S{Q, k). A relevant subset of admissible meshes can 
be dehned according to the result of Proposition 9 in [6] . 

Definition 6. Let Q he the mesh of active elements defined according to o and 
([2]) with respect to the domain hierarchy D D ^ . A mesh Q is 

strictly admissible of class m if 


cu 


i—m+l 


(3) 


where 


^l-m+l _ y I g . g ^ ^ ^ ^ + 1) C } , 

for £ = m,m + l,...,N — 1. 

The overlay Q* of two meshes Qi, Q 2 is the mesh obtained as the coarsest 
common rehnement of Qi and Q 2 , usually indicated as 

Q* = Qi® Q2- 

Let {flf}£=o,...,Afi-i and {f22}^=o,...,7V2-i be the nested sequence of domains that 
dehne the hierarchical meshes Qi and Q 2 , respectively, with = ^ 2 - The domain 
hierarchy }£=o,...,v,-i, with = max(A^i, ^” 2 ), associated to Q* satishes 

U ^2 and (uf 3 wf U 

for £ = 1,..., — 1, where fif = 0 when i > Ni, for i=l,2. Consequently, if Qi 

and Q 2 are strictly admissible, for any level £, we have 


since any Q & either belongs to or and the overlay Q* is a rehnement 
of both Qi and Q 2 . The overlay Q* of two strictly admissible meshes is then a 
strictly admissible mesh. Note that the number of elements of the overlay mesh 
Q* is bounded as follows, see e.g., [H |22], 

#Q. = #(Qi ® Q 2 ) < #Qi + #Q 2 - Qo, 

where Qq is the initial mesh conhguration. The above inequality may be used 
for discussing the optimality of the resulting adaptive isogeometric method, anal¬ 
ogously to adaptive hnite element setting. 
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2.3 The REFINE module 


In order to define an automatic strategy to steer the adaptive method, we will 
propagate the rehnement procedure over a suitable neighborhood of any marked 
element. 

Definition 7. The neighborhood of Q E Q (1 with respect to m is defined as 
M{Q, Q, m) :={Q'e :3Q"e S(Q, £ - m + 2), Q" C Q'j , 

when £ — m+1 >0, and M{Q, Q,m) = % for £ — m + 1 < Q. 

A sequence of strictly admissible meshes can be recursively dehned by suitably 
extending the rehnement of coarser regions beyond the set of marked elements M 
through the algorithms presented in Figured! 

Note that these algorithms follow the structure of informal high-level descrip¬ 
tions in the spirit of the analogous modules related to the adaptive hnite element 
methods. 


Q = REFINE_RECURSIVE(Q, Q, m) 

for all Q' € Af{Q, Q, m) do 

Q = REFINE_RECURSIVE(Q, Q', m) 

end for 

subdivide Q and 

update Q by replacing Q with its children 

Figure 1: The REFINE and REFINE_RECURSIVE modules. 

By exploiting key properties of the REFINE_RECURSIVE module, summa¬ 
rized in Lemma 0 and Proposition [9] below. Corollary [10] characterizes the output 
of the REFINE procedure [6]. 

Lemma 8 (Recursive rehnement). Let Q be a strictly admissible mesh of class m 
and Q E Q. The call to Q* = REFINE_RECURSIVE(Q, Q, m) terminates and 
returns a refined mesh Q* with elements that either were already active in Q or 
are obtained by single refinement of an element of Q. 

In addition, ii Q E the level £* of all newly created elements Q* generated 
by the call to Q* = REFINE_RECURSIVE(Q, Q,m) satishes 

£* <£+l. (4) 

In order to verify this, we may note that the recursion is applied to elements of level 
< £, and, in particular, of level < £ — m + 1. If Q* is a child of Q then £* = £ -f 1. 


Q* = REFINE(Q,7W,m) 

for all Q G Qn A4 do 

Q = REFINE.RECURSIVE(Q, Q, m) 

end for 

Q* = Q 
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Otherwise, Q* is obtained by splitting some elements in the chain of neighborhoods 
generated by set of recursive calls and, consequently, i*<i — m + 2<i + l since 
m > 2. 

Proposition 9. Let Q be a strictly admissible mesh of class m > 2 and let 
Q G , for some Q < i < N — 1. Then it follows that the call to Q* = 
REFINE_RECURSIVE(Q, Q, m) returns a strictly admissible mesh Q* ^ Q of 
class m. 

Corollary 10. Let Q be a strictly admissible mesh of class m >2 and M the set 
of elements of Q marked for refinement. The call to Q* = REFINE {Q,Ai,m) 
terminates and returns a strictly admissible mesh Q* ^ Q of class m. 

Note that in each computation of the neighborhood Q, m), the choice of 

level i — m + 2 for the support extension yields the smallest neighborhood that is 
acceptable for preserving the class of admissibility of the mesh when subdividing 
the given element Q. Nevertheless, depending on the underlying hierarchical mesh 
conhgurations, the basis functions could be also truncated at different intermediate 
levels. 


3 Linear Complexity 

This section is devoted to a complexity estimate in the style proposed by Binev, 
Dahmen and DeVore [3] and, in an alternative version, by Stevenson [23], for 
adaptive hnite element methods. 

3.1 Auxilliary results 

For every pair of mesh elements {Q, Q'), let dist((5, Q') be the Euclidean distance 
of their midpoints. Given a Q & all Q' G ffiQ, Q, m) satisfy 

dist((5, Q') < diam(S'((5, £ — m + 2)), 


where I = I{Q) and 

diam(^(Q, I-m + 2))-.= 2-^+™-2 (2p+l) = 2-^C ,, 
with Cg = Cs{p,m) := 2™“^(2p+ 1), p := ^Pi- Hence, 

dist(g,g') < 2-^-^Cd, Cd = Cd{d,p,m) := ^/dCg. (5) 
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Lemma 11. Let Q be a strictly admissible mesh of class m > 2, M. the set 
of elements of Q marked for refinement, and Q' G Q fl Al. Any newly created 
<5 G Q*\Q obtained by the call to Q* = REFINE_RECURSIVE(Q, Q',m) satisfies 

dist{Q,Q') <2-^^^'>C with C := VdC, (7:= (2-^ + --(6) 

1 _ 2^ j 

where then C depends on d,p and m. 

Proof. The existence of Q G Q* \ Q means that REFINE_RECURSIVE is called 
over a sequence of elements Q' = Qj,Qj-i,...,Qo and corresponding meshes 
Qj, • • •, Qo so that Qj-i G Af{Qj, Qj, m), with Q' E Ai and Q being a child of Qo, 
namely £{Q) = £{Qo) + 1- Since £{Qj-i) = £{Qj) — m + 1, h follows 

£iQj) = £{Qo)+j{m-l). (7) 


We have 

dist((5,Q') < dist((5,(5o) + dist((5o,Q0 

and 


./ 

dist((5,(5o) = 2~^^'^h~^Vd, dist(Qo,Q') < dist(Qj, Qj-i) . 

i=i 

According to (jS]) and ([7]), we obtain 


5^dist(Q„Q,_i) < ^ 

j=i j=i j=i 

°° o-^iQo) 9-qQ)+i 

- Ch. = z -- C,.. 

j=0 


< 


1 - 21 - 


1-2 


1—m 


Hence, dist(Q,(5') < 2 dQ)c^ where C is the constant dehned in (jb]). □ 


3.2 Main result 

The main result of this paper states the existence of a generic constant A = 
A{d,p,m) < oo that bounds, for any sequence of successive rehnements, the ratio 
between the number of new elements in the final mesh Qj and the number of all 
marked elements encountered in the rehnement process from Qo fo Qj. 





Theorem 12 (Complexity of REFINE). Let M. := U/=o of marked 
elements used to generate the sequence of strictly admissible meshes Qo, Qi, • • •, Qj 
starting from Qo = G^, namely 

Qj = REFmE{Qj_i,Mj-i,m), /or j e {1,..., J} . 

Then, there exists a constant A > 0 so that 


j-i 

#Qj-#Qo< 

j=0 

with A = A{d,p,m) := 4(4(7 + l)"^, where C is defined in (j6|). 

Proof. We denote by G := [J^- the set of the initial mesh elements and all 
elements that can be generated from their successive dyadic subdivision. Let 
Q G G, Qj\,4 G AA, let 

r 2TQ)-mM) if + 1 and dist(g, Qm) < C, 

Qm) •= \ 

[0 otherwise. 

The proof consists of two main steps devoted to identify 

(i) a lower bound for the sum of the A function as Qm varies in Ai so that each 
g G Qj \ Qo satishes 

A(g,gA,) >1; (8) 

Qm^M 

(ii) an upper bound for the sum of the A function as the rehned element Q varies 
in Qj \ Qo so that, for any j = 0,...,J—1, each Qm G Aij satishes 

A(g,g7n) <A. (9) 

QgQj\Qo 

If inequalities ([H]) and ([9]) hold for a certain constant A, we have 

#Qj-#Qo = E IS E E -^(g, Qm) 

Q&Qj\Qo Q&Qj\Qo Qm^M 

, 7-1 

< 5^ a = a5^#ai,, 

Q J —0 

and the proof of the theorem is complete. We detail below the analysis of (i) and 
(ii). 
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(i) Let Q G Qj \ Qo be an element generated in the refinement process from 
Qo to Qj, and let ji < J he the index so that Q G Qj^+i \ Qji- Lemma fTTl together 
with dl]) state the existence of Qi G Aij^ with 

dist(Q,Qi) < and £{Q) < i{Qi) + I, 

and, conseqnently X{Q,Qi) = > g. The repeated use of Lemma [m 

yields a sequence {Q 2 , Qs, • • • } with Qi_i G Qj^+i \ Qj^, for ji > j 2 > js > ■ ■ ■, 
and Qi G Aij^ such that 

dist((5i-i, Qi) < C and i{Qi-i) < £{Qi) + 1. (10) 

We iteratively apply Lemma [H] as long as 

X{Q, Qi) > 0 and £{Qi) > 0 , 

until we reach the first index L with X{Q,Ql) = 0 or £(Ql) = 0. By considering 
the three possible cases below, inequality ([H]) may be derived as follows. 

• If £{Ql) = 0 and X{Q, QQ > 0, then 

51 ^(Q. Qm) > A(Q, Ql) = > 1 , 

since £{Q) > £{Ql) = 0. 

• If X{Q,Ql) = 0 because £{Q) > £{Ql) + 1, then flTU]) yields £{Ql-i) < 
£{Ql) + 1 < £{Q) and hence 

X{Q,Qm) > X{Q,Ql-i) = > 1. 

• If X{Q,Ql) = 0 because dist((5, Ql) > C, then a triangle inequality 

combined with Lemma [m leads to 

L-l L-l 

2 i-hQ) c < dist(Q, Ql) + ^ dist(Q„ Qi+i) < C+ J2 C . 

i=l i=l 

Consequently, < YldZi ^^^1 we obtain 

1 < = 5;)a(q,o.) < 55 kq,qm)- 

i=l i=l Qm&M 
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(ii) Inequality ([9]) can be derived as follows. For any 0 < j < J—1, we consider 
the set of elements of level j whose distance from Qm is less than C dehned 
as 

According to the dehnition of A, the set B{Q_M,j) collects the elements at level j 
so that X{Q, Qm) > 0- We then have 

hQxi+i 

^ ^ A(0,Qa<)= Y1 (11) 

QgQj\So QgG\So 

Since the diagonal of an element Q of level j is 2~^ \/d, the diagonal of the 
hypercube composed by the union of the closure of all elements in B{QmG) is less 
or equal to 

2 • 2^-^ C + 2-^Vd = 2-^ \/d (4 C + 1), 
where C is dehned in ([H]). Hence, 

(4(7+1)", 

and the index substitution k ■.= 1 — j + (.{Qm) reduces flTTll to 

Q^QAQo ^=0 

oo 

< 2J2‘^-'^#B{Qm,j) = di^B{QM,j) < A, 

k =0 

with A = A{d,p, m) = 4(4(7 + 1)". □ 

4 Conclusions 

We developed a complexity estimate which says that the ratio between the rehned 
elements and the marked elements along the rehnement history stays bounded, 
when rehnement is performed as proposed in [6]. In particular, this estimate 
guarantees that the rehnement routine ensuring the (strict) admissibility of the 
rehned mesh remains local at least when looking at the overall rehnement process. 
Note that for a single rehnement step, it may be impossible to prove such an 
estimate [23] . 

Our work paves the way to the analysis of optimal convergence of the adaptive 
strategy proposed in [B] that will be addressed in further studies [7]. 
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